Signatures of molecular correlations in the few-electron dynamics of coupled quantum 

dots 
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We study the effect of Coulomb interaction on the few-electron dynamics in coupled semiconductor 
quantum dots by exact diagonalization of the few-body Hamiltonian. The oscillation of carriers is 
strongly affected by the number of confined electrons and by the strength of the interdot correlations. 
Single-frequency oscillations are found for either uncorrelated or highly correlated states, while 
mult i- frequency oscillations take place in the intermediate regime. Moreover, Coulomb interaction 
renders few-particle oscillations sensitive to perturbations in spatial directions other than that of the 
tunneling, contrary to the single-particle case. The inclusion of acoustic phonon scattering does not 
modify the carrier dynamics substantially at short times, but can damp oscillation modes selectively 
at long times. 
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Low-dimensional heterostructures enable direct prob- 
ing of the time evolution of carriers. In particular, charge 
oscillations between coupled quantum structures have 
been measured in a number of different systems, such as 
Josephson junctions quantum wells Q, or quantum 
dots [3j, |4(. Coupled quantum dots (CQDs), where the 
number of confined carriers can be controlled experimen- 
tally, are a most interesting case as they posses a discrete 
energy spectrum which stems from the quantum confine- 
ment in all three spatial directions, thus constituting the 
physical realization of the particle-in-the-box problem. 
Understanding charge oscillations in these structures is 
not only of fundamental importance but also of techno- 
logical relevance, since control of the coherent dynamics 
of charge states is at the base of many proposals for novel 
nanoelectronic devices [H, [f| 0] and quantum logic gates 
[1> S E3] • This has motivated recent attempts to achieve 
control of charge localization and dynamics in CQD sys- 
tems by means of microwave excitations [ll|, static Q and 
time-dependent 12, HI electric fields. 
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Coherent charge oscillations between CQDs have been 
recently demonstrated in the single-electron 11[ and few- 
electron [3, IH regimes. The origin and behavior of these 
oscillations in the former case is well understood: when a 
single electron is placed in one of the dots of a CQD sys- 
tem, with no other electron in it, it oscillates back and 
forth between the dots with a frequency given by the 
energy difference between the bonding and antibonding 
"molecular" states. When a larger number of electrons 
is present, the system behavior is much less understood. 
Recent experiments working in the latter regime exposed 
a single- frequency oscillation of the carriers [H, 0| , which 
was interpreted in terms of an effective single-electron 
picture. However, in general one would expect more com- 
plicated oscillation patterns owing to the non-trival den- 
sity of states of Coulomb-correlated few-body systems 



In this work we theoretically investigate charge oscilla- 
tions in few-electron CQDs. We show that, as an effect of 
Coulomb correlations, the amplitude, period, and shape 
of these oscillations are strongly dependent on the num- 
ber of electrons confined in the structure. Either single- 
frequency or multi-frequency charge oscillations occur, 
depending on the strength of the correlations between 
molecular levels of the CQD. Additional external mag- 
netic fields are shown to provide a versatile means for 
tuning the period of the charge oscillations. This unique 
behavior due to the full spatial quantization is character- 
istic of quantum dots and clearly differs from the plasmon 
oscillations in coupled quantum wells [1, 17]. We also in- 
vestigate the effect of the acoustic phonon bath, and show 
that the different strength of electron-phonon interaction 
for different few-electron states leads to selective suppres- 
sion of modes in the multi-frequency dynamics. 

We consider a system of two vertically coupled dots, 
as depicted in the first panel of Fig. [TJ populated 
with N interacting electrons. The dots are gated 
GaAs/AlGaAs heterostructures, as those built in Ref.ll8, 
The conduction-band single-electron states are described 
within a three-dimensional envelope function approxima- 
tion, including an electric field E z along the z direction 
19]. For this cylindrically symmetric configuration, the 
single-particle eigenfunctions can be given the separable 
form Vw ffs (r, er) = 4>nm(x,y)<p g (z)Xa(o), with n and m 
radial and angular quantum numbers of the Fock-Darwin 
state [l6j], g — 0, 1 labeling the bonding and antibonding 
eigenfunctions associated with the double-well potential, 
and s the spin orientation. In order to compute the few- 
electron states exactly, we use a full configuration inter- 
action approach [13] : the single-particle states are popu- 
lated with N electrons in all possible ways to construct a 
basis of Slater determinants \§>j), where j stands for the 
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set of many-body quantum numbers. Then, the three- 
dimensional TV-electron Hamiltonian is diagonalized. We 
first assume a closed system, so that coherent charge os- 
cillations take place. In the last part of this work, we 
briefly investigate the effect of phonons. 

Initially, the CQD is subject to an electrostatic bias E z 
which favors localization in the lower dot (QD2). Then, 
at time t = 0, the bias is removed (non-adiabatically) and 
the carriers start oscillating between the two dots. In or- 
der to simulate this process, the configuration interaction 
calculation described above is performed twice: first by 
taking E z > 0, then E z = 0. The computed eigenstates 
are = £.,• cy|$j), and |* ; ) = £j c yl $ j>> with and 
without the E z field, respectively. We take, as the ini- 
tial state of our few-particle system, the ground state in 
the biased condition, |Wo). Its time evolution, after the 
removal of the bias, can be computed from 

l*o(*)> = 5>-* B,t (*j|*o)|*J>. (1) 
l 

Here the sum runs over the I unbiased states, whose en- 
ergy is Ei, participating in the spectral decomposition 
of |^o)- 111 order to visualize the charge oscillations, we 
evaluate the particle density inside the upper dot (QDI), 
Ni. This magnitude has been resolved experimentally 
using pump-and-probe techniques Q or single-electron 
transistors [J]. 

In Fig. Q] we illustrate the time-dependent occupancy of 
QDI for N = 1, 2, 3, 4, 5. Gray and black lines represent 
weak and strong initialization bias (see caption), respec- 
tively, as schematically reported by the two V(z) profiles 
in the first panel. For N = 1, we retrieve the expected 
sinusoidal oscillation, whose frequency is given by the 
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FIG. 1: Schematic of the CQD structure under study and 
occupation probability in the upper dot as a function of time 
for iV-electron systems. Gray and black lines represent the 
oscillation patterns for weak (E z = 50 kV/m) and strong 
(E z = 170 kV/m) initialization biases, respectively. The dots 
have thickness W = 10 nm, lateral harmonic confinement 
Tiojq = 3 meV and the interdot separation is Lb — 8 nm. The 
prominent oscillation period depends mainly on the number 
of oscillating electrons (note the two different scales for the 
time). 




FIG. 2: Occupation probability in the upper dot as a function 
of time for iV-electron systems. Gray and black lines represent 
the oscillation patterns with (B = 5 T) and without (5 = 
T) magnetic field, respectively. The CQD structure and biases 
E z are the same as in Fig. [Tj 



bonding-antibonding energy separation, A^bab- Both 
weak and strong biases completely localize the electron in 
QD2, thus giving rise to essentially the same curve. How- 
ever, as N increases, the initialization bias starts playing 
a critical role, as it determines the number of electrons 
localized in each dot at time zero. This, in turn, gives 
rise to very different oscillation patterns (compare black 
and gray curves for N > 1). 

An inspection of the few-electron oscillations in Fig. [T] 
shows that they exhibit multiple frequencies with differ- 
ent amplitudes (most apparent in the N = 2 panel). This 
is a signature of Coulomb interaction in the interdot dy- 
namics, which can be understood from the coefficients of 
the spectral decomposition of the biased state, (^i\^ ): 
only the unbiased states which have a finite overlap will 
contribute to the time evolution in Eq. ff]). In the single- 
electron case the in-plane and vertical degrees of freedom 
decouple. The initial state is localized in QD2 and results 
to be the linear combination of one bonding and one an- 
tibonding unbiased states. This is because the in-plane 
xy component of the wave function is the same with and 
without bias, while, in z direction, the ground biased 
state is given by ipo{z) +y?i(z). Therefore, only these two 
states contribute to the spectral decomposition and the 
oscillation frequency is given by their energy difference. 
By contrast, in the few-electron case Coulomb interac- 
tion mixes the radial and vertical degrees of freedom, so 
that the in-plane parts of the wave functions with and 
without E z are no longer identical. As a result, the spec- 
tral decomposition may involve several pairs of bonding- 
antibonding states with different radial wave functions. 
Each pair contributes to the time evolution with its own 
frequency, given by its bonding-antibonding energy split- 
ting, and an amplitude that is proportional to its spectral 
decomposition weight. 

The mixing of radial and vertical degrees of freedom 
due to Coulomb interaction renders charge oscillations 
sensitive to perturbations in the xy plane, even though 
they do not affect the single-particle tunneling. This 
is shown in Fig. ® where we compare single- and few- 
electron oscillations in the presence and absence of a 
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FIG. 3: (a) Main unbiased states contributing to the spectral 
decomposition of the ground state in the N = 2 CQD of Fig.Q] 
at E z = 150 kV/m. Occupation probability of the upper dot 
in systems with: (b) null (e* = oo), (c) partially quenched 
(e* = 10 3 ), and (d) regular (e* = 12.9) Coulomb interaction. 



magnetic field applied along z. While the N = 1 os- 
cillation is insensitive to the field, the frequencies of the 
N = 2 and N — 3 ones are increased. This signature 
of electron-electron interaction could be used in experi- 
ments to distinguish between the dynamics of indepen- 
dent and correlated electrons. It also offers a unique way 
to modulate the few-particle oscillation frequency. 

To gain further insight into the role of Coulomb inter- 
action, in Fig. [3] we compare the oscillation patterns in a 
CQD with no interaction (panel b), artificially quenched 
interaction (panel c), and regular interaction (panel d). 
In all cases, the electric field initializes both electrons in 
QD2, and yet major changes take place as the Coulomb 
intensity is varied. Again, the changes can be understood 
by analyzing the main unbiased states contributing to 
the spectral decomposition of the biased state. These 
are the s-shell singlets, ^ and ^3, represented in 
Fig. El^a) by their dominating electronic configuration. 
In the schematic representation, electrons are indicated 
by arrows, whose direction is the spin state, and can oc- 
cupy the single-particle bonding (lower horizontal line) 
or anti-bonding (upper horizontal line) states. For van- 
ishing Coulomb interaction, ^Sq is composed with equal 
weight of the two bonding-antibonding pairs + ^2) 
and (^2 + ^3)- The frequency arising from the two pairs 
is identical, &E™ B = AE* 3 AB = E t (E t is the single- 
particle tunneling energy), and the oscillation shown in 
Fig-E^b) is reminiscent of a single-particle one. Switching 
on a weak Coulomb perturbation [Fig. \3[c)] introduces a 
small departure from this limit, here A£^ B = 0.24 meV 
and A£^ B — 0.20 meV: now the two pairs have simi- 
lar amplitudes but different energies. The superposition 
of the two harmonic motions gives rise to a beat with 
modulated amplitude and frequency. Finally, for a real- 
istic Coulomb interaction [Fig.[3](d)], A£^ B = 1.26 meV 
and A£^ B = 0.04 meV, the weights of the two bonding- 
antibonding pairs are very different. As a consequence, 
there is one dominating oscillation mode with large am- 
plitude and low frequency, and another with small am- 
plitude and high frequency. 




FIG. 4: Occupation probability in the upper dot as a function 
of time for N = 2 (upper row) and N = 3 (lower row). Here 
E z = 500 kV/m and W — 10 nm. The lateral confinement 
and interdot barrier thickness in each column are: (a) hujQ — 1 
meV and Lt = 6 nm; (b) htoo = 6 meV and Lt = 6 nm; (c) 
Tiujo = 1 meV and Lb = 10 nm. e* = 12.9 in all cases. 



The sizable changes of the bonding-antibonding en- 
ergy splittings in the presence of Coulomb interaction 
are an effect of the electronic correlation between molec- 
ular states, often disregarded in previous studies of 
multi-particle dynamics in coupled quantum dots 13j, [14 1 
and wells [l3]- The weaker the interdot correlation the 
stronger the multi-frequency character, and vice-versa. 
Therefore, one can control the nature of the charge os- 
cillations by designing CQD structures in the regimes of 
either weak or strong interdot correlation with respect to 
Et ■ This is shown in Fig. [5J where we compare the N = 2 
and N = 3 charge oscillations for different structural pa- 
rameters, keeping a realistic value for the Coulomb inter- 
action throughout. In (a), the interdot barrier is thin and 
hence the tunneling energy is large. The large splitting 
between molecular orbitals implies weak interdot corre- 
lation, which leads to an oscillation pattern resembling 
the multi- frequency beat of Fig. [3jc). In (b), the radial 
confinement is increased, which increases the Coulomb 
repulsion within the structure. Since the vertical con- 
finement is unchanged, interdot correlation moves into 
the strong regime and a quasi-single-frequency behavior, 
similar to that of Fig.[3jd), is retrieved. In (c), the barrier 
is made thicker with respect to (a) , so that the tunneling 
energy diminishes. This again enhances the interdot cor- 
relation, leading to the quasi-single-frequency behavior. 
We point out that the small tunneling energy could be 
the reason for the single-frequency oscillations reported 
in Ref. [3J, U experiments. 

The electron dynamics in CQDs is severely affected 
by dissipative processes [2l[. In what follows, we shall 
investigate which novel effects may appear in the few- 
electron charge oscillations due to the interaction with 
the 20 mK acoustic (deformation-potential and piezo- 
electric) phonon bath. We calculate the transition rate 
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FIG. 5: Occupation probability as a function of time in a 
N — 2 and TV = 3 CQD subject to phonon interaction. The 
insets zoom-in at short and long times. Note that the TV = 2 
case evolves into a single-frequency oscillation. The CQDs 
have TujJq, W = 10 nm, Lb = 6 nm (TV = 2) and Lb — 8 nm 
(TV = 3), and are initialized with a bias E z = 500 kV/m. 



7ii between the correlated states and \^>i) accord- 
ing to Ref. and employ the Pauli master equa- 
tion to evaluate the system's time evolution: = 

j-{Ej - Ei)pij - J2i llj V H Pi] + % 12i7uPU- We stress 
that here pij is an element of the density matrix corre- 
sponding to the (unbiased) few-particle states, with en- 
ergies Ei and Ej. 

Figure \E\ shows the resulting charge oscillations for a 
CQD with N — 2 (left panel) and TV = 3 (right panel). 
In both cases the oscillation amplitude is clearly damped 
by the phonons. However, since the transition rate is 
different for each couple of initial and final states, the 
different modes contributing to the multi-frequency oscil- 
lations are damped at different rates. As a consequence, 
in some CQD structures all but one mode are quickly 
suppressed, and the initially multi-frequency oscillation 
turns into a single-frequency one at later times (compare 
the insets in the TV = 2 panel for short and long t). For 
TV > 2, a higher number of correlated states participate 
in the spectral decomposition of the biased state, so that 
more oscillation modes show up in the charge oscillation. 
As a result, it is difficult to find conditions where only one 
mode survives the phonon damping (see insets of TV = 3 
panel). Therefore, this is unlikely to be responsible for 
the single-frequency oscillations observed in Refs. 0. 

In summary, we have shown that the few-electron dy- 
namics of CQDs is strongly affected by interdot electronic 
correlations. Drastic changes in the oscillation pattern, 
from simple sinusoidal to complicated beats, take place 
depending on the correlation strength. Electronic inter- 
action further renders few-electron dynamics very sensi- 
tive to perturbations in all directions of the space, and 
not only in that of the tunneling. Upon inclusion of 
phonon damping, few-electron charge oscillations may 
experience an additional effect, namely a selective sup- 
pression of frequency modes. This is however unlikely 
to be responsible for the single-particle-like oscillations 
reported in early experiments 0, 0| , which can be under- 
stood as an effect of strong molecular correlation, due to 
the small tunneling energy. 
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